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SUMMARY 


Th.is report describes and explains the computer program that 
has been used to solve the governing equations of the potential 
flow in the cross-sectional planes of a radial inflow turbine 
scroll. The derivation of the governing equations and the 
description of the numerical solution can be found in NASA CR 
(to be assigned) Report entitled "Analysis of the Cross Flow in 
a Radial Inflow Turbine Scroll." 

A list of the main program, the subroutines, and typical 
output example are included. 
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COMPUTER PROGRAM 


The Fortran computer program, which is used to solve the 
equations of potential flow in the cross-sectional plane of 
the centrifugal machine's volute or scroll will be de- 
scribed and listed. The input data to this program includes 
information about inlet flow properties, the cross section 
geometry, and its corresponding boimdary conditions, which are 
shown in Figures 1 and 2. The maximum allowable number of mesh 
lines in the x and y directions are 50 and 30, respectively. 

The program output includes an echo print of the input data. 

/ 

At each grid point, the value of the potential function, the 
velocity components in the x and y directions, as well as the 
magnitude of the velocity vector and the angle it makes with the 
X axis are printed. Additional output is obtained in the form 
of plots- of 'Constant potential function contours, and arrows 
showing the Idirection of the velocity in the axial-radial plane. 

The logical relations of the main program and siibroutines 
are shown in Figures A-1, A-2 and A-3. A brief description of 
the main program and siabroutines is given below. 

Main Program ; 

The main program is primarily used to call the subroutines . 

The logic flow diagram is illustrated in Figure A-1. 

Subroutine IRL (J) ; 

It is called from the main program and the siibroutines to 
determine IL(J) and IR(J) , the first and last interior mesh points 
for the scroll cross-section at a given J mesh line (see Fig. A-4) . 

Subroutine INTPL(I) ; 

It is called from the main program and the subroutines to 
determine JL(I) and JU(I) , the lov7er and the upper interior mesh 
points for the scroll cross-section ah a given I mesh line (see Fig. A-4) . 

Subroutine MESH (I, J) ; 

It is called from the main program and the subroutines to 
identify the interior mesh point neighbors, calculate the coeffi- 
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cients of the five points Laplace difference operator, and to 
determine the most recent values of the potential functions at 
the interior mesh points . 

Silbroutine SHABMF ; 

It is called from the main program, to determine the source 
strength at every mesh point F(I,J) , for a given mass flow rate 
and through flow velocity profile. The logic flow diagram is 
illustrated in Figure A-2. 

Subroutine AWATEF ; 

It is called from the main program to determine the values 
of the potential function, at each grid point. The logic flow 
diagram of this subroutine is illustrated in Figure A-3. 

Subroutine BCSI(I) ; 

It is called from the subroutine AWATEF, to determine the 

values of the potential function at the boundary points on a 

given I mesh line. . ORIGINAti PAGE IS 

OF POOR QUALITY 

Subroutine BCSJ(J) : 

It is called from the subroutine AWATEF, to determine the 
values of the potential function at the boundary points on a 
given J mesh line. 

Subroutine HWEL ; 

It is called from the main program to determine the values of 
the velocity components in the x and y directions, at every mesh 
point. 

Subroutine ABDLAH ; 

It is called from the main program to calculate the magnitude 
of the velocity vector in the axial radial plane and the angle it 
makes with the x-axis . 

Subroutine HOSHY ; 

It is called from the subroutine SHABAN to calculate the 
flow density at each grid point. 

Sxibroutine PLTISO ; 

It is called from the main program to plot the contours of 
constanr potential function. 
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Additional library subroutines namely, PLOT (XI, X2 , X3) , 
SYMBOL (Yl, Y2, Y3 , Y4 , Y5 , Y6) , and ENDPLT, are called from 
the main program, and used in the plotting of both the potential 
contours and the velocity direction. 
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FIGURE (A~2) SUBROUTINE SHABAN 













DEFINITION OF SYMBOLS AND FORTRAN LISTING 


Symbols 

B 

DO 

EPSMAX 


F 


GAMA 

H 

ICOM 


IFG, IPAP, 
IFREVO 

IP 


IPLT 


ISYM 


Description 

Width of the scroll exit (Figs 1, 2) . 

Scroll inlet diameter. 

The largest tolerable value of the square of the 
sum of the absolute values of the deviations of 
j from their previously computed values. 

Vector containing the values of Poisson’s source 
strength at each grid point. 

Ratio of specific heats. 

Length of scroll exit nozzle. Figs. 1 and 2. 

Parameter to control type of calculations . 

ICOM = 0 for incompressible flow and ICOM = 1 
for compressible flow. 

Parameters to control the mass source distribution 
(i.e., the type of through flow velocity profile) - 

Frequence of intermediate printouts. — Solutions 
are printed after every IP iteration. 

Nximber of contours in the output plot for the 
potential velocity. 

Parameter to control type of scroll cross section. 
For symmetric cross-section ISYM = 1, and for 
unsymmetric cross-section ISYM = 0. 


ITLIAX Maximum allowable number of iterations . 

IVPLOT Parameter to control plotting of program output, 

IVPLOT = 1, a plot is prepared; if IVPLOT = 0, 
no plot IS prepared. 


Ml, M2, M3, N3 


Number of mesh lines of the scroll cross-section. 
Pig. (A-4) . 


ORIGINAL 
OF POOR QUAUP^i 


ORIGINAL PAGE IS 
OF POdfe QUALinf 
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Symbols 


Description 


P 

PBX^PBY 

P^(i=0,l,2,3, 

R 

RC 

Rl 

RG 

ROA 

ROBX,ROBY 

THETA 

TO 

U,V 

DBX,UBY 

VBX^VBY 

VE 


Vector containing the values of the potential 
function at all grid points i 

Vectors containing the values of <j> on the scroll 
boundaries . 

The values of the potential function at the standard 
five points Laplace difference operator. 

Radius of the scroll cross-section equal D/2, 

Figs . 1 and 2 . 

Radius of exit portion of the unsymmetric cross- 
section, Fig. 3. (Ref. 1) 

Radius of scroll exit - 

Gas constant. 

Vector containing the values of the density at all 
grid points . 

Vectors containing the values of the density at 
the boundary points. 

Vector containing the values of the angle betv/een 
the velocity vector V, and the x-axis at each 
grid point. 

Stagnation temperature 

Vectors containing the values of the flow velocity, 
components in the x and y directions, respectively. 

Vectors containing the values of the velocity 
component in the x- direction at the boundary 
points of all I and J mesh lines. 

Vectors containing the values of the velocity 
component in the y-direction at the boundary points 
of all I and J mesh lines. 

Scroll exit velocity in radial direction. 
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Symbols 

VISOBR 

W 

WBX.WBY 


WO, WOP 


Desoription 

Numerical values of the velocity potential con- 
tours to be plotted as output. 

Vector containing the absolute values of the 
velocity vector V, at each grid point. 

Vectors containing the absolute values of the 
velocity vector V, at the boundary points for 
each I and J mesh lines, respectively. 

Successive relaxation factors. WO for interior 
mesh points, and WOP for the points adjacent to 
the boundaries . 
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PROGRAM LISTING 


IV G I.HVEL 21 MAIN DATt = f7ZS}3 21/43/41 _ _ _ 0001 

” ’ ■ IMpJlCIT ^EAL*e (A~H,0-Z> 

DI MANSION X(5 0,21) ,1T(50,21 I . VAR( 50 • 21 ) , VRA130. 21 ) 

CCMMON /ONE/ VI 

COMMOV /TaO/ ll_0 t JV.R , IbVM, IPAP. 1 cod. ip. 1 com 

COMMON /=OJR/ Rl.RF.IFHEVO 

COMflOV /DIMS/ R. RC ,d .OX.DY.Bl ,3B 

CCMWO.V /FDRM/ Ml . M2, IC . JC. II C, JJC, M3.N3 

C CKMON /MS>^/ S 1 ,32 ,33 ,S4 ,P1 ,P2 .P3.^4 , JL, JU, IL, Ih 

COMMON /Pril/ =>{ 50, 301 ,PdYC 2,30) ,P3X( 50,2) ,PPC5 0,30 J ,?3Y1 12 ,30) , 

1 PSXl (50 ,2) , FCsO.SO ) 

COMMON /VHES/U (5 0 , 3 0) , V{ 50 , 30 ) ,UBX{ 50 ,2) ,VBX(50 .2) , 

I Do Y (2, 30) , V3YC 2,50 ) ,U1( 50 ,30) , VOt 5 0,30) .03X1 (50,2) 

1 . Vaxl C50,21,UBY1(2,30),V3YH2.30),VV(50,30) ,ANol_E (5 0, 30) 

CCHMJV /CAi-/ M21 , I Cl , JCl , Ml 1 ,N31 . I C2 
common /3H/ RG.GAMa , TO.MV, VH 

common /mo/ EPSMAX ,EPSMX2, WC.MOP.ir max, ITMX2 
COMMON /R3/ HOA( 50 ,30) ,R03X( 50,2) , ROBY (2.30) 

COMMOM /go/ VV SX ( 50, 2),VVBY (2,30), ANGl. KS0.2) , ANGU2( 2, 30 J 
CALL UND=LW 
C 

CC CALCULATE 2-01 M . P3TEN510NA- FUNCTION G VELOCITY VECTOR V IN THE 
CC SCROLL CROSS SECTION FOR INCQMP. G COMP. FLO* 

C 

CC READ INPUT DATA #*:4**4TT*4*****4*«**S:»**<^4^*»4A»t*T**444;*»-^;<i*.» 

C 

REA0(5.102) RG .3AMA, T&.VE 
R£AD(5.I00) Ml ,M2. DO ,R1.3.H,R,RC 
RE AD (5,1 01 ) ICGM.ISYM, IFG , IPAP.IFREVO 
REA0(5,200) ITMAX,EP5MAX,ViO,»QP 
R£AD(5.7177) I P. IV P^OT . I P_ T 


ITMX2= UMAX 
EPSMX2==PSMAX 


Vl=l . ODO 
RC=RC/( 2<=R > 

E=B/(2*R) 

IF( I 5YK.EQ. 1 ) B = B/2.0D0 
h=H/(2»-I) 

?0=R1/D0 

R1=R1/(2TR) 

* V=16. OJO#R 0?34R/D0 
IFdSYH.EQ. I ) *.V = 2.000TtoV 
RF=2.000*!R 
R=R/(2*R) 

RATIC=(R1+H+R)/(2.000*R) 
DX=2*R/( Ml- 1 ) 

DY=2FR/( M2- I ) 
IFUSYH.EQ.l) DY=R/(M2-1) 
a I=0X/3Y 
BB=81*3I 
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IV G L:V=L 21 


MAIN 


DATt = 77233 


Zl/Ao/Ml 


M2l=H2-J 

1 C-={ Ml +1 )/2 

JC=(M2-M)/2 

r F ( I 3YM .Ea. 1 ) JC=1 

N3 = ( D/'i)Y-0. OOOIOD3 )+i 

V3=( (S+a+^J /DX-O. 0 00 1000) n 

ICI= IC-1 

JC1=JC-1 

1C2=1C+J 

N31=N3-1 

M3l=M3-l 

IFClSYM.hQt 1) GO TO 555 
I I c=t»<c/ox~o.oouiooo)+Mi 
J JC^ US+RO/OY-O.OOOJ 0D0) + i 

c 

555 00 12 1=1. M3 
DO 12 J=l.M2 
angle { 1 . J)=0 00 
RCAC 1, JJ=1. OOO 
ROBXt 1,11=1 .000 
rtOBX (1 .21=1 .000 
ftOBYl 1 , J)=1.000 
ROBY (2. J 1=1 .000 
F II , J)=D.0D0 

12 continue 

Ir(lFG.EQ.O) GO TO 5557 

00 3552 I=2,M31 
CALL INTPLI I ) 

5552 READ (5, 5553 ) (FU ,J 1, J=Jw JUJ 
5557 C0NTINU5 
C 

CC » RITE J>J = JT DATA FOR CriEC< = = 

c 

WRITE! 6,1313) 

WkJTc(6,405) ^G, GAHA.TQ, V£ 

WRI TE I 6,34 3 ) Ml , M2 ,Kl ,B.H . RC . RF, DO 
WRITE! 6.310 ) ICOM, 1SYM.1F3, IPAPjIFRcVO 
wRITE!6,320 ) I TM AX . EP5MA X, i’D. WOP 
WRITE!6,33C,> I P, I VPLOT.IPLT 
WKITc!6,30i) OX.DY.IC, JC,M3,N3 
C 

CC SOLVE FOR INCOMP. FLU» 

C 

1*R=0 
I COD=0 
CALL SHA3AN 
CALL A«ATEF 

1 LD=0 

CALL -IVVEL 
CALL AB0-A1 

IFUCOH.EO.O) GO TO 313 


PAoE U00«. 


original n® K 
OF POOR 
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N IV G level 21 MA _ DAJE_= 77£S2 . - - 21/A.S74-I _ _ - - PAI7E’ Ot) 03 

CC SULVE = OR COM®. FL3W **¥*♦#** + + 

c 

DO 2 1=2. M3 
CALL I NT=Lt 1 ) 

DO 2 J=JL.JU 
2 PPCI.J)=1.0D0 

VEC=VE 
17=0 
I n f?=I 

717 IT=IT+1 

ERRaR=O.ODO 

IC0D=1 

V£=VcO/'H03Y 12,2) 

CALL SHA3AN 
CALL AXATEF 
ILO=0 

CALL rtVVEL 
call A30_AH 
DO 7 I =2. M3 
call INT®L(I) 

DO 7 J=JL,JO 
SAV=PP(I ,J) 

7 E«RQR=ERRaR+OABS (® ( I > J )-SA V J »-=!:2 

DO 111 1=2. M3 
CALL INTPlII) 

30 111 J=JL,JU 

111 P®C .J ) = ®U ,J) 

DO 112 I =2, M3 
30 1 12 J=l. 2 

112 P3X1 £I , J J = PBXJ I, J) 

00 114 J=2.K21 

30 114 1=1,2 
114 pan <1. J) = ®SY( 1, J) 

IF{IT.G£.1THX2 ) GO TO 91? 

IF{£RRDR.GE-EPSMX2) GU TO 717 
919 (.'?UT£(6. 14) 
feRITC(6,407) 

DO 9 1=2. M3 

9 HrPlTEIo, 408) 1,(H( I, J).J=I,M21) 

WRITEC6.151) IT, ERROR, VE 
WRITE(6,A04) 

DO 40 J=1,M21 

40 1»R1TE(6,400) J ,( P( I , J) ,1=1 ,M30 

1F{ ISYM.H3. 1) GO TO 46 
VkRITE<6.431 > 

0 0 4 2 1=2,K3 

42 WRITE! 6, 402) I , { P3X{ 1 , J) . J = 1 , 2 } 

GO TO 45 

46 WRITE! 6, 4005) 

DO 47 1=2, M3 

47 WRITE!o,406) I,PBXI1,2) 

45 WRITE!6,403) 

DO A4 J=l,M2i 
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IV G LHVhL 21 MAIN DATs = 77283 21/45/Al PAGE 

a* *KITE<6,a02) J , ( P3''{ i . J ) . I = 1 . 2 ) 

I »P=0 

CALL riVVEL 
CALL A30_AH 

8ia continue 

IF n VPL3T.EQ.0 ) GD TO 7007 
C 

CC “LOTTING 0= THE FLO* OIRECTIDNS *^=^:*A:*:*A:<=»^***W* **:t***!tti»»! 

c 

ISKIP=1 
JSMI P=1 
0 AD=0. 25 
CALL SEG=_T 
ICPM=2 
ICP=1 

00 59 1=1, M3 
0 0 5S J=1 ,M2 
XU . J)=( I-n *=0X0 
Y ( I, J) = { J-I I^DXO 

58 CONTINUE 

59 CONTINUE 
V.KITei6, la J 
WRITE! 6,5001 5 
DC 62 1=1, M3 

62 ,RITE<6,2C) ( X ( i , J J , J=l, M2 J 

teRITECO.lA) 

WRITE! 6,5002) 

DO oa 1=1, M3 

6A WRITEI6.20) (Ytl ,J) , ,M2) 

WRITE! 6, lA) 

75 IF ! I SYM.EQ. I ) GO TO 76 
call plot I 1 . 5, I- 5,23 ) 

GO TO 77 

76 CALL PLOT U .5, 3.3,23 ) 

77 I0P0DW=1 
SI2E=0.1 A 
I SYMBL=75 

IFilCP.EO. 1 ) 30 TJ II 
IF!ICP.E0.2) GO TD 133 
II DO 50 J=2,H3 
DO 50 J=1,M21 
VRK! I , J )=270.0 +ANGi_E ( I, J ) 

50 CONTINUE 

CALL PL3TIO.0,O.0.5J 
00 501 1=2, M31 
CALL INTPLII) 

DO 511 J=JL,JU 
XXC=XII . J) 

YYC=Y!1 , J) 

CALL SYMBOL { XXC , Y YC , SI ZE . 1 S . V RR ( I , J ) ,-l ) 

511 CONTINUE 
5C1 CONTINUE 

CALL P_3T !0. 0,0. 0.3) 


®ginal page is 

©f POOR QUALITY 
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IV G LEVEL 21 MAIN DATE = 772o3 ai/AS/Al- - PAGE 'O&O'S ‘ 

I CP= ICP+1 

IF(I CP.GT . ICPM ) GO TO 3 
CALL PAUSE 
CALL PAUSE 
GC TO 75 
133 CONTINUE 
C 

CC CONTOURS Pl_OT + sis* 

C 

DO 510 1=1, H3 
DO 510 J=1 . M2 
510 VARI I . J)=P{ I ,JJ 
DO 51A 1=2, M3 
CALL INTPLII) 

J1=JL-1 

J2=JU+I 

VAR{ I , J2} = PBX( 1,2) 

IF(I SYM.EQ.l) GO TO 514 
VAR< I, J1 ) = ?BX< 1, 1 ) 

SIA CONTINUE 

00 515 J=2,M21 
CALL IR^CJ) 

1 1 =I L- 1 
I 2=1 H+ 1 

V API 11 , J ) = =3Y( 1, J ) 

515 VARI 12 , J) = PBY(2. J) 

CA_L PLOT 10.0,0.0,3) 

DO 222 I=_=l,IPuT 
READ 15, 5) VISD5R 

CALL PLT ISO! M2 1, 1 S YMBL , I U=OOV. , SI ZE , V I SOBR . VAR , X, Y) 

CALL PLOT (0.0.0.0,23) 

222 CONTINUE 

CALL PLDT 10.0,0.0,3) 

3 CONTINUE 

CALL END=LT 
7007 CONTINUE 
5 FORMAT IF 13 .0 ) 

14 FORMAT I 1-il) 

20 FORMATIIOI 6X,=6.3) ) 

ICO FCRMATI2I 10 ,6F10 .3 ) 

101 F0RMATI3I10) 

102 FORMAT! 4F10 .0) 

151 FORM ATI 23 X , • 1 T CO.M= • , I 5 , 7 , 20 A, *eRR02= > . E2C ,S,/,2CX, *V EX. = *,E20.&) 

2 00 FORMATtI10,3F10, Oj 

301 FORMAT! 20X, ‘OX = • ,F 10. S , 5X , • O Y = • ,F10, 5 , 5X , ’ I C =‘,110./. 

1 20X.*JC =’,113, 5X,‘M3 = • . I 1 0. 5X, • V3 =’,110,/) 

310 FORMAT!/, 20X, * ICUM = • ,i 1 0 , / , 2oX . • 1 5 YM =’, 1 I 0 #/ ,20X, • IF G =*. 

1 I 10, /, 20X, » I = AP =» ,I10,/,20X,«IFREVO=’ ,I 10,/) 

320 FORMAT! /. 2 OX, • ITMAX=' ,I10,/.20X. ’EPS MAX=» , FI 0.8, /. 20X, ’»0 =• , 

1 F 10 .5, /, 20X ,* V.OP =*,= 10.5./) 

330 FORM AT {/ , 20X, • IP = • , 1 1 0 . / . 20X , • I VPlCT= • , 1 1C ./ ,20X, • I PLT =• . 

1 110 ,//) 

340 FORMAT! /,20X, • Nil = • , I 1 0 , / ,20X, *M2 =' , II 0 , / ,20X , • R1 =',F10.S 
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5.N I V 


G L.5VEI. £1 


KA IN 


DATE = 77263 




I t/.20X.*B =• . = lO.S./i20X.'K =• ,F10.5./.20X,‘ rfC =•. 

1 FlC . Jt SOX , • O = • ,r I 0 ,5.2Xt • JNl T DIM. ' , 7 ,20A.‘ JO =«.FIC.b 

1 .2X,'UMT DIM.’ ./) 

400 FORMAT!//. 7X, • J = • . I 2. /, 1 0 1 5X, 501 5. 6/ J I 

401 FORMAT (//. 15X I ’ . SX, • P 3X 1 • . laX. ’P8X2 • J 

4 02 FORMAT!/, I OX. I 4, 5X ,2 ( 5 X.Dl 5 .6 ) J 

4C3 F0R.^AT! //, 1 5X, *J ‘ ,5X, 'PaTl * . I5X, 'PBY2 •> 

4t-4 FORMAT!//, 20X, *VA_JES OF Vt^DCITY POTENTIAL FD-i COMPnESSIbO 

l£ FLOW’./J 

AOS FORMAT (//.20X, *GAS CONSTANT =• . = 10. 5,/,20X, ‘GAMA =*,FlC-3 

1 ,/.20X,« stag. TEMP. =• ,rlO .5 ,/,20X, ’RAD. EXIT VE_= > , F 1 0 .5 , /) 

406 FORMAT! /, 15X, I 4, lOX.O 15.61 

407 FORMAT !/.l OX.' PD ISSONS SOURCE 0 1 S. AT EACH GRID POINT',//) 

406 FORMAT! //, 7 X.' X =• , I 1 0 ,/ .4 ( SF 1 O .3 1 1 

1313 FORMAT! //, I OX, ’IN3UT DATA',//) 

400S FORMAT !// ,20X, «I • . 5X, • PBX2 • . / ) 

500t =ORMAT!//,20X.»X-CCDRO. OF THE PLOTTlNo DOMAIN',//) 

5002 FORMAT!//, 20X. 'Y-C30R0, 0= THE PLOTTING DOMAIN',//) 

7177 FORMAT !3I 10) 


5553 F0RMAT!3!F10 . J )) 
7177 FafiPAT!3I10) 

STOP 

END 


PA^E 0006 


ORIGINAL PAGE IS 
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IV G LiVEL 21 A1HATE.F - DATE = 7 72E3 21/A5/A1 PAGc 0001 

StlBRDUTINE AWATEF 
IMPLICIT REALMS (A-H»D-Z) 

common /one/ VI 

COMMON /T»0/ ILD , I teR. ISYM. IPtP . ICOJ t IP»ICOM 
COMMON /DIMS/ R. RC» M» B, DA »Oy »B1 . S3 
CQM'JON /FORM/ Ml .M2 • IC.JC. I IC. JJC«M3 *N3 
CCMMOM /MSH/ si. S2tS5t54.=>l tP2.P3.=>A» JL.JU.i:-, 1R 

COMMON /PHI / P(SO»30> ,PBY (2.30) . PdX < 50.2).PPC50.30).P3yil2.3C). 

J P3XK 50, 2J.r( 50.30) 

COMMON /CAL/ M21 . ICl . 3C1, M 1 1.N31. I C2 
COMMON /HO/ EPSMAX,CPSMX2.*O.»0P.iTMAX,lTMX2 
C 

C CALCULATE POTENTIAL VELOCITY 3 1 ST . IN THE SCHOLL CROSS SECTION BY 
CC SOLVING POISSONS EON. Vt I TH' NEU MANN O. C. THE ITER. METHOD IS S.O.R. 

C 

2 OO 1 1=1, M3 
DO I J=1 .M2 
1 Ptl,J>=0,000 

ITER=0 

15 1TER=ITER+1 
EPS=0.000 
OO 5 1=2, M3 

xxi=iic-n »OX 

Xl=OASS(XXl ) 

YY1=DA3S (R*R-X 1=X1 ) 

Y1=DSQHT( YYl) 
call. INTP_(I) 

DC A 0=JL,JU 
AL2W0 

Wl=l .ODO-W 

*lB=W/ (2.DDC+2,ODO*BB ) 

00=1 . 000-W3*^B3 
wai=l .ODO-UfS 
YY2={ JC-J)*DY 
Y2=DABS( YY2) 

XX2=0A3S( R*R-Y2* Y2 ) 

X2=DSQRTtXX2) 

SAVP=P(I ,J) 

CALL IR_(J) 

CALL HESHU.J) 

LI 1=1. ODD 

L33=1.0D0 

WBI1=1 .000-»3<=L1 1 

»S 12=1 . 000- i*3*L33 

C2A=1/{0Y*0Y*S2=SA) 

C13=1/(DX«'0XAS1* S3) 

S2A= 1/t DY aDY*( SH+SA) ) 

SI3=l/IOX*OX*(SltS3) ) 

S 1L=513*L1 1, 

S1K=S13*L33 ■' 

Pll=Lll*Pl/S’l I 

P22=P2/S2 

P23=L33*P3/S3 
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P4A=P4/S4 

IF( I .cQ.Ml . and ,J .EO.JU ) GD TO 923 
i F ( t • ^ 0 • ll.* AND-* J • JE } GO TO BL> 

IF<I .£Q.1_.AND.J .EQ.JU} G3 TO 83 
IFtS.EQ.IL) GO TO 90 
IF(I .£0,1?.AN0.J.EQ. JL) GO TO 91 
IFil .EQ.lrl.ANO.J .Sj.JU) G3 TO 92 
IFU .SQ.IR) GO TO 93 
IF{ J .EO .JL) GO TO 94 
IF(J.EO.JU) GD TO 9S 

P (1 . J)-=( 1 . 000- WOP) *P I I .J) F{V.CJP/ (2 + 2*86 ) )«(L1 1*P 1+L33*33 + 

1 aQwtP2+=4)-ra iJ)*DX«=DXl 
GO TO 477 
8C CONTINUE 

1F(J.£Q.1) GO TO bl 
Pai=Y2»Bl*P( I. J+l )/X2 
Pg2=Xl*^( I+l ,J)/(3l*YI ) 

CKl-{ W*S1_» ( ( l/Sl)-Y2*Sl/X2) +W*S24*( ( I /S4j -XI/ (B1*Y1) ) )/lC13+C24) 
IF(CK1 .GT .5 ,999000 .AND. C< I .uT . 1.001 ) GO TO 16 
CK=1.0D0-£;X1 
GO TO 17 

16 CXI=CKl/w 
CN=1,000-1K1 
VH = O.ODO 

W=3 . OOO 

17 CONTINUE 

P(1,J)=(;<1»P(1,J) + a1c(S134C =31+P33) +S24*( P22+PB2)-F < 1 . J) 72) 

I /IC13+C24) )/CK 
GD TO 477 
31 CONTINUE 

P<l,J)j=(WlcP(I ,J > + W6» 12*BB*P2 + P33-F( 1 , J)«DX40X J l/WEll 
GO TO 477 
85 CONTINUE 

P31=Y2*31»P{ 1 , J-n 7X2 
pe2=xi*p( 1 + 1 . j )7(3 i*Yi ; 

CK1 = (W-31L*< ( 1 731 )-Y2»B17X2 ) + W*524*( { l/S2)-X 17 (Bl»Vl) ) )/(C13+C24J 
1=(CK1. 3T. 0,995000, AND. CKl »UT.l. 001 J GO TO IS 
CK=1 ,0D0-CKI 
GO TO 19 

18 CX1=CKII/V, 

CK=1 .ODO-CKl 
tel=^O.ODO 

W= 1.000 

19 CONTINUE 

PC I , J) = ( iH*P( I ,J)+W* tS13*( P8l+P33)+3 24*{P44+PB2)-FCl, J }72) 

1 7{ CI3+C24) }7CK 

GO TO 477 
90 CONTINUE 

IFCJ.GT.JC) GO TO 90S 
PB1=PCI.J+1)*Y2#317X2 

CKl=CVf*SlL*( 1 17Sl)-Y24bl 7X2) )/(C13+C24) 

IFiCKl ,3T .0 .99SODO .ANO,C<l.L-T. 1. 001 ) GO TO 20 
CX=1 .ODO-CKi 
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GO TO 21 

20 CKl=CKl/h 
CS:=1 .ODO-CK 1 
V> 1=0 .COO 

W= 1.000 

21 CONTINUE 

P(I.JI=(#1*P(I,J )+W4 iS13«{PBl +P33)+S24»(P22+P44J-FU . J) /2) 

1 /{ C13«-C24) )/CK 

GO TO 477 

909 CONTINUE 

P81=Pl I. J-1 1*Y24D1/X2 

CKl=l4*51L»< { I/al J-y2*Bl/X2) )/t C13FC24 ) 

IFCCKl .GT.0.9990DD. AND. CKl .LT.l. 0011 GO TO 22 
CK=1 .OOO-CKl 
GO TO 23 

22 CX1=CK1/H 
CK=1 .ODO-CK 1 
>*1=0.3 00 
>*=1. ODO 

23 CONTINUE 

P(l. JJ-={»1*P(I .J J + «N>lS13*tPai+P33H'324-*( P22+P44J-F (I . JJ/2J 
1 /1C13+C24) J/CS. 

GO TO 477 
51 CONTINUE 

IFIJ.EQ.l) GO TU 910 

CK1 = I tV.*SlK/53) + ft*S24/S4) / {C13+C24) 

IF (CKl . GT. 0.599003 .AND. CKl .1_T .1 .001 1 GO TO 24 
CK=1 .OOO-CKl 
GO TO 25 

24 CKl=CKl/»/ 

CK=1 .ODO-CK 1 
W1=0 .003 
W=l. OOO 

25 CONTINUE 

P{I»J) = {i*l*P(I,JJ + !if« (S13*P1 1+S1K4VI *DX+S24*?22-F a .J) / 2J 
1 /(C13*C24) J/C<. 

GO TO 477 

910 CKl=W*Sl</tS3*(C13+C24)> 

IFICKI.GT. 0.999000. AND.CX1.lt. 1.001) GO TO 26 
CK=l .OOO-CXi 
60 TO 27 

26 CK1=CK1/Vi 
CK=1 .303-CKl 
1*1 = 0.000 
rf=l .000 

27 CONTINUE 

P{I,J)=(>*1*P(1 ,J)*-W4=IS13*P1 1+S1K*V14DX+2*S24*P22-F U. J)/2) 

1 /{C13+C24) )/CK 

GO TO 477 
92 CONTINUE 

IFtJ.GT.V3> GO TO 921 

CX1 = ( C W*S1K/S3 ) + H*S24/S2)/ t C13 + C24) 

IFICXl .GT.O .99S000.AND.CK1 .lt. 1. 001 ) GO TO 28 
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ci<;=i . ooo-:k 1 
GD TO 29 
2E CKi=CKl/* 

CK=i .odo-:k 1 
»I=C .200 
K=1 .030 

29 continue 

P<I«J) = (*l*PlI.J)4-i»i*{5I3*?ll + SI< »VlTi)x+S24»P4A-P 11 fJ)/2J 
1 /IC13+C24) )/CX. 

GO TO 477 

921 continue 
IFUSYM.EQ. n GU TO 922 
JFtJ.GT.JJC) GO TO 922 
XCl=RC-t 1-Ml )*0X 
YYC1=0A3S( RC4RC- XC14XC1 3 
YCI=DSQRT{YYC1 3 
YC2=(RC+81-( J-1 ) *0Y 
XXC=DA8S( 4C*RC-YC2*YC2) 

XC2=DS0RT(XXC) 

PS2=XC1 *P( £-1 , J1 /IS14YC13 
= B3=YC2*Sl*P (1 ,J-13/'XC2 

CK1 = <W*3 244{ li /3 2)-XCi/< B1*YC1 3 )+WSSlK*£ £ IXSB 3-YC243 1/XC 23 3 /( C 13 

1 

IFICKI .GT .0 .999000. AND. C< 1 .CT. 1 . 001) GO TO 30 
CK=1 .3D0-CK1 
GO TO 31 

30 CKl=C<l/» 

CX=I . ODO-CKl 
w 1=0,000 
W=1.0D0 

31 CONTINUE 

P{ I . J)=£ »1*P{ 1 , S134 IP83+P1 13 I-S24* £ P44 + PB2 3-F £ I . J 3/2 3 

1 /( C13-I-C243 3/C< 

GO TO 477 

922 CONTINUE 
PB3=Y2*B1*P£ 1. J-13/X2 
PB2=Xl*^P£!-l .J)/£314Y13 

CK1 = (W*SIK4£ £ l/S33-Y2*Bl/X23 + «»’S24T£ £ 1 /S2 )-Xl / £ 3 1 * Y1 3 3 3 /£Clo + C24 ) 
3F£CK1 .Gr ,0,9993DD,AND.C<1 .wT, 1 .001 ) 60 TO 32 
CK=1.000-C<1 
GO TO 33 

32 CKI=CK1/H 
CK=1,ODO-CK1 
k 1 = C .000 

»=1 .000 

33 CONTINUE 

P£l.J3=(nl*P(I.J3 + »*(313»( PS3+P1 13 + S24»tPA4-t-PB23-F£I .J3/23 
1 /(CI3 + C243 3/C< 

GO TO 477 
93 CONTINUE 

IF£J.GT.N33 GO TO 930 
CKl=w*Sl</£ S3* £C13+C243 3 

lF(Ci<l .GT. 0.999003. ANO.C<l.l_T.l. 001) GO TO 34 
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CK=1.000-C<1 
GO TO 35 

34 CKl=OCJ/w 
CK=1 .003-C<1 
u 1 = 0.000 

w = 1 ,000 

35 continue 

P(I,J>={Wl*PtI ,J)«-W*IS13*P1 l+S 1N*V1 *DX+S24*(P22+P44J-F ( 1 > J 1/21 
1 ✓{ Cl 3+C24 1 )/C< 

GO TO 477 

930 CONTINUE 

IP(ISY«. = O.H GO TO 932 
IF(J.GT.JJC) GC TO 931 
yC2=(HC+3)- ( J- 1} «OY 
XXC=OA3S( RC*RC-YC24YC2> 

XC2=D5QRT (XXC) 

PS3=YC2«Bl*PtI .J-:)/XC2 

CKI={W4S1X*U 17S3)-YC2*B1/XC2) )7<C13+C24) 

IPICKl .GT .0 .9990 03 ,AND.C< 1 .-T .1 .00 1 1 GO TO 3o 
CK=1 .000-CKl 
GO TO 37 

36 CK1=CK1/* 

CK=1.000-EK1 

a1=0,C09 

» = 1 . 000 

37 CONTINUE 

P(I,J>-=I*1AP( I.J )+W»{ S13*( P33+P1 1} <-S24T( P22+P44J-F (1 , J)/2J 
1 /( C13+C24 ) J/CK 

GQ TO 477 

931 CONTINUE 

IFIJ.GT. JCIGC TO 932 

= Cl.J> = I*tl«P(I.J }4-»3*( PI U-33AI P22+P44)-FI I . J1*0X4UXJ )/HB12 
GO TO 477 

932 CONTINUE 
PB3=Y2*81^P( I. J- 1)7X2 

CKl = IV.*St K4{ (I/S3)-Y2=»Bl/X2))/(C13+C24) 

IF (CXI. GT. 0.999000. AND, CXI .LT. 1.001 ) GOTO 3S 
CX=l .ODO-CX 1 
GC TO 39 
3£ CKt-=CXl7W 

CK=1 .300-C<1 
K 1 = 0.000 
»= 1 . 03 0 
39 CONTINUE 

P(I.J) = tWl*P(I,J)FW*(SI3*( P33 + P! I JFS24*(P22 + P4 4)-F( 1 , J J/2) 

1 7t C13*C24) )/C< 

GO TO 477 
94 CONTINUE 

IFIJ.EQ.IJ 60 TO 945 
IF(I.GE.IC) GQ TO 940 
P34=x!i*P(H-l,JJ/(al4YI) 

CK1 = («I*S24*( tl/5 4)-Xl/(31*Yl) ) )/(Cl3FC24) 

IF (CXI .GT .9 .9990 03 .AND. CO ._T .1 .001 ) GOTO 40 
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CK=1 .DDO-CKl 
GO TO 41 

40 CK1=C<1/* 

CK=1 .0D3-CK1 
w! = v^.ODO 

1 .000 

41 CONTINUE 

»(I.J)=(W1*PII,J11W4(S154C?33 + P1 m-S24*{P22+PB4 )-F II, J 1/21 
1 7(C13+C24) )^C< 

GC TO 477 
94S CONTINUE 

P<l,J)=rfl+P(I.JH-*S*< PI 1 +3 33 + 2. 0D0768*P22-F( I ,3)*UX43X1 
GO TO 477 
940 CONTINUE 

P<I,J} = l»I4PU,J)+*3*( Pll+333+BB»?22~F( I .J J*0X*DXJ )/33 
GC TO 477 
95 CONTINUE 

IFd.GE.lO GO TO 950 
PS2=X1*P(I+1 .J)/(3I4Y1 1 

CKi = < W*S244U ly'S2J-XI/lbl4 YD ) 1 / (C13+C24J 
IF (CKl .GT.S .999300 .AND. C<1 ,LT . 1 ,001 1 GO TO 42 
CK=1 .ODO-CKl 
GO TO 43 

42 CK1 = CK1('» 

CK=1 .ODO-CKl 

K 1=0 .000 

fc=l .ODO 
45 CONTINUE 

Pll,Jl=(414P(l,J)+w*(S13*( P33+P1 ll+S24*{P44+PB2J-Ftl ,J)/23 
1 /(CI3+C24) }/Cv 

GO TO 477 
925 CONTINUE 

P{ I .31=0.000 
GO TO 477 

950 CONTINUE 

1F(I .GT.MI ) GO TO 951 
P32=X1*PU- I ,J)/131#Y1 1 

CM=tW*S24*({ I/S2)-XI/I3I*Yin )7<C13+C24J 
IFICKl .GT. 0.999000 .ANO.CKl ,LT . I .001 ) GOTO 44 
CK=1 .ODO-CKl 
GO TO 4b 

44 CK1=CK1/H 
CK=1 .ODO-CKl 
Vi 1=0 .0 03 
W=l. 000 

45 CONTINUE 

P{I.J) = <vili«^P(I,J ) + »#=(S134lP33+PI 1 )+S244CP44+PB2)-F (1 . 3 1/2} 
1 /t C13+C24) l/CK 

60 TO 477 

951 CONTINUE 
IFII5YM.EO.il GO ID 952 
IF(I.GT.IIC) GO T3 952 
XC1=RC-U-M1 )*DX 
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A*.AT = ® 


rVCr=DABS'(RC4RC-XCl*=XCl ) 

YC:=D50aT(VYCI ) 

Pe2=XCl*3{ I-I.J}/<31*YC1J 

CKI=rf¥S24*t n/S2)-XCl/(fal*YCl ) )/<C13+C24) 

IF ICKl.&T. 0.9950a0.AND.CiU .l-T. I .001 ) GO TO 4o 
CK=1 .000-CKl 
GO TO 47 
Ato CKlxCKl/w 

CK=J .0O3-C< 1 
hl=0.000 
1 .000 

47 CONTINUE 

P( I . J) = (K1*P(I lS13*lPll+P33)+324*IPAA+PB2J-F{i. J J72) 

1 ✓( C13+C24) ) /CK 

GO TO 477 
952 CONTINUE 

CKl = W«S24/< S2A(C13+C24n 

IF (CKl .ST .C .9993DC .ANO.C< I .-T .1 .001 } GO TO 43 
CK=1.0D0-CK1 
GO TO 49 
4S CKl=CKl/* 

CK=1.0D0-CK1 
h 1=0.000 
W=l. 000 
49 CONTINUE 

PCI.J)=l*l*PtI»J )+»»( bl34( P II+PSBJ+SEA^PAA-FII .JJ/2) 

1 /ICI3+C24))/CK 

477 CONTINUE 

ePS=EPS+DA3S(0( 1 ,J )-5AVP)=4^2 

4 CONTINUE 

5 CONTINUE 

1FI1TER7IP4IP .NE. ITERJ GO TO 2222 
WRITEI6.102) EPS, ITER 
DO 900 J=1.M2 

90C hRITEIS. 10 1 ) (PU. J 1. l=l.vl3 ) 

WRITE! 6. 1043 
00 1000 1=1, M3 

lOC" WRITEt6,I03) I ,(PBX!I,J),J=1.2) 

WRITE<6, 1053 
DO 1130 J=l.M2 

1100 »RlTe(6,l03) J.{P3YI1,J3.1=I,2) 

2222 CONTINUE 

:FIE°S.wT.1 ,0003 go to 3 

WO=O.90DO«WC 

WOP= 0.9030* WO® 

GO TO 2 

3 IFCITER .GE. ITMAX) GO TO 3 
IFIEPS .GE. EPSMAX) GO TO 15 
s continue 

WRITE <6,1023 EPS. ITER 
WRITE! 6. 143 

DC 6 I =2. M3 
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CALL I NTPL( 1 t 

CALL acsn I ) 

6 CONTINUE 

iDQ 7 J=1 , M2 1 
call IRL(J) 

CALi_ SCSJ ( J. VI ) 

7 CONTINUE 
IFtlwR.EO.U GO TO 222 
WRITE(tj.l06 ) 

DO 9 J=I ,M2 I 

9 tiRITEIo, lOl ) J .( P( I , JJ ,1 = 1 ,M3J 
I F( ISYM.EG. 1 ) GO TO 333 
V.R1TE(6, lOA) 

DO 10 1=2, M3 

1C *RIT£(6,1Q3) I .{PBXCI.JJ ,J=1,2J 
GO TO 

333 ttRITEIo, 107) 

DC 12 I =2, M3 

12 »RITE(6,103) 1,PBXII.2J 

AAA ARITE16.135) 

0011 3=1 ,M21 

11 1*R1TE( 6, 103) J ,(P3Y{ I ,J) .1 = 1,21 

1 A FORMAT (/// ■ • *♦ »*==»==A**-=**A=A*«**=A***»A**4iA**A=*= 

1 *♦»*«♦♦***»** ♦ A *=«*==**» = *»!:**»♦*♦** **=*#=*=» *m**^****^**Vt:*lf ‘ ,// J 

101 FORMAT! //,7X, • J = • , 12,/, 10! 3X,3D1S.6/) ; 

1C2 FORMAT!//, lOX.'EPa = • , D1 2 . 4/ 1 0 X. • NO . OF ITERAT IONS=* , 1 A) 

1 03 FORMAT!/, 1 OX.IA, 5X,2!bX,015.5}» 

lOA FORM AT (//, 1 SX, • I • ,5X, 'P3X1* ,1 5X. 'PBXa' J 

ICS FORMAT!//, 1 SX, *J • ,5X, • PBYl • ,1SX, • PSVa • 1 

106 FORMAT! //,20X, *VA_U£S OF POTENTIAL VELOCITY P!X,J) FOR INCOMPHESS 
HOLE F_0**./1 

107 FGRMATt//,20X, *I ' ,5X , ■ P3X2* ,/ J 
103 FORMAT!/, 15X.1 A, 10X.D15.6) 

222 CONTINUE 
RETURN 
END 
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SU3RDJTINE A3D_Ari 
IMPLICIT HEAL*6 {A-ri,Q-2) 
common /ONE/" Vl 

COMMON /'T»3j' IUD . IkR • 1SYM» IPAP, IcOJ t IP»1CUM 
common /form/ mi ,M 2, iCt JC,I IC, JJC. M3, N3 
CUM'oON /MSH/ S 1* S2 ,53 ,S4 .PI .P2.P3.PA . JL, JO, IL. IF 
COMMON /VEES/O<50, 30 J ,Vl 53.3a ).UBX( aO , 2) . V8X{ 30 .2) . 

I US Y( 2. 301. V3Y( 2. 30 1,01(50 >301 , VO (50 >30 1 .UbXl (50,2) 

1 , VBXU 30.2) .UBYK 2.30) . VBY 1(2.301 ,Vy(50.301 .ANGLE (50 .30) 

COMMON /CAL/ M2I . ICl • JCl > Ml 1.N31.I C2 

common /so/ VVBX(30 ,2 1 .V VBY12 .30) . ANGLl (50.21 . ANGL2( 2, 30 1 
C 

CC CALCULATE THE VELOCITY VECTOR V IN THE SCHOLL CROSS SECTION BY 
CC UjING THE TWO COMPONENTS U.V IN THE X&V DIHECTIOWS 
C 

IF(IWR.EO.l) GOTO 5 
*RITE(6.14) 

WRITE! 6. 149) 

WRITfclS.lSO ) 

5 CONTINUE 

OD 2 1=2, M3 
CALL INTPL(I) 

DO 2 J=JL,JU 
U2=DABS( UII . J) J 
V2=DASS(V( I, J) ) 

VV(I ♦J)=OSORT{ L2=U2+V2*V2) 

IF(U( I. J l.GT.O .ODO.AND .V( 1 , J) .LT.O, ODO) GO TO 82 
IF(U( 1 . J) .LT .0 .0 03 . AND.V ( 1 . J ) .GT .0 .000 ) GO TO 53 
IF(U( I ,J) .LT.O. ODO.AND. V(I . J) .LT.O. QOO) GOTO 84 
1F(J { 1 . J ) .£0.0 .ODO.AND.V( 1 . J) .GT.O. 000) GO TO 85 
I F (U ( 1 . J) .£0.0 .003 .ANO.V ( I , J ) .LT .0 .000 ) GO TO 86 
IF(U(t .J). ST. 0.000. AND. V( 1 .J) .EQ. 0.000 ) GO TO 87 
I F(U (1 > J ) ..T.O .000 .AND. V( 1 .J) .£0.0.000 ) GO TO 88 
1F(U( 1 . J) .EQ.O .000 .AND ,V (i . J) .EG.O. 000) GO TO 37 
ANGLE! I.J) = (DATAN( V2/U2) ) *57 . 296000 
GO TO 1 

82 ANGL£(I.J)=(6.233000-(DATAN(V2/U2}) )*57.29oODO 

GO TO I 

S3 ANGLE! I, J)= (3. 1423 DO- ( DAT AM! V2/U2) 1 ) *57.296030 
GO TO 1 

84 ANGLE( I.J)=(3. 142000+ (OATAN! V2/U2) >)*S7.296000 
GO TO 1 

65 ANCLE! I . J)=90. 000 
GO TO 1 

86 ANGLE! 1 ,J)=270 .303 
GO TO 1 

57 ANO.EC I.J >=0.000 
GO TO 1 'I 

88 ANGLECI . J)=180.0D0 

1 CONTINUE ‘ 

IFdWR.EQ.l ) GO TO 2 

kRIT£!6.100) I.J.VVtI.J), ANGLE(I.J) 

2 CONTINUE 
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WkITE(6.200» 

DC 3 1=2. M3 
DC 3 J=l,2 
J2=D4SS( JBX( I. J) ) 

V2=DABS (VBX ( I . J) ) 

VVEXt 1 . Jl =DSQST( U2*U2+V2*V2J 

IF {D3x ( I , J ) .GT . 0 .OOO.AND. V3X( 1 , J} .-T. C« 000} GO TO 72 
IFtjSXd .J) .LT .0 .ODU .AND. VBX( 1 .J }. ST .0.000) GO TO 73 
IFIuaXd . J) .LT.O.0D0.AND.V5XI I , J) .-T.O.ODO J GO TO 74 
IFtUSXt I, J ) .E3 .O.ODO.ANO.VBXI I .J ).GT. O.QDO) GO TO 75 
IF{U3X( I , J) .EO.O .OJO . ANU.VBX ( I , Ji.UT .0.000) GO TO 7o 
lF(uaX{ I . J J ,GT .0.030. AND, VBX( I . J) .£0. 0.000) GO TO 77 
IFIUSXI I. J ) .EQ .0 .:03 .ANO.VBXC 1 ,3 ).£3. 0.000) GO TO 77 
ANGLI a .J1=(0ATAN{ V2/U2))*57. 296000 

1F{U3X{ I, J) .LT .O.ODO.ANG.VBXU . Jl.EQ. 0.000) GO TO 76 
GO TO 13 

72 ANG1.1 (I . J) = (3. 23 30 00— (OAT AN (V2/02) ) )«57. 296006 
60 TO 10 

73 ANGLIU . J)= (3. 1423 00- tOATAN(V2/L»2) ) ) *57.296030 
30 TO 10 

74 ANGLI ( I . J )= {3, 1420O0+COATANI V2/02J ) ) *57.296030 
GO TO 13 

75 AN6L II 1 . J) = 90. ODO 
GO TO 10 

76 ANGLI II .J)=270 ,000 
GO TO 10 

77 ANG^ 1 I I .3 )=0.0 00 
GO TO 10 

78 ANGl 1(1 ,0 = 130.000 
10 CONTINUE 

IFIIlkR.EO.l ) GC TO 3 

aRITEI 6. 100> I .J . VV3XI I .J) .ANGLI (1 .3) 

3 CONTINUE 

WRtTE<6»300) 

DO 4 J=2.M21 
DO 4 1=1,2 
J2=DABS(UBy(I. J) ) 

V2=UABS( V37I I . J) ) 

VVBY ( I . J ) = OSQR T{ U2*J2+V2*V2) 

IFIUBYI 1 .J) .GT .0.300. AND. VBYII .Jl.LT .0. ODO) GO TO 62 
IFIUBYI I > J) .LT .0.030.AN3.V5YU , J) .3T. C.ODO) GO TO 63 
IF (U3YI 1 , J) .LT .6 .309 .ANO.VSYd.J ).-T .0.000) GO TO 64 
IF (UBYd ,J) .EU. 0.333 .AND. VBYd .J) .GT. 0.000) GO TO 65 
IFIUBYI I. J) .£0 .0.000. AND. V3Yd ,J).LT. 0. ODO) GO TO 6b 
IFtuBYd ,J> .GT .0 ,306 .AND. VBYd.u ),ea . C.ODO) GO TO 67 
IFIUBYd .J) .LT. 0.000. AND. VBYd .3>.cQ.O. ODO) GO TO 68 
IF luSYI I. J ) .EO .0.000. AND. VSYI 1 .JU.EQ. 0.000) GO TO 67 
ANGL2 I I, J >= IDA TANIV2/U2)) *67.296003 
GO TO 20 

62 ANGL2 (I. J)= I6.2630D3-I DATANl V2/U2) ) ) *57-296000 
GO TO 20 

63 ANGL2I 1 , J)=I 3. 1420DO-IDATANI V2/u2) ) )* 57.296000 
GO TO 23 
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bA 'anGU 2{ I , J } = (3. l-»23O0 +< DAFAnIC V£yU2J } } *57.29600 0 
GC TO 20 

65 ANGL2C I . J)=9C. COO 
GC TO 20 

60 ANGL2CI , J)=270.000 
GO TO 20 

67 ANGG2( I . J}=0 .0 DO 
GO TO 20 

66 ANGL2( I. J)=ISO.OOO 
2C CONTINUE 

IFUHfH.EQ. X> GO TO A 

V«FITE(6.100) I .J« VV3V( 1. J) . ANGL2II . J ) 

A continue 

I A =CRNAT<1-H) 

100 FOPVAT I20X. IS. SX» IS. SX.E20.9. SX.FS.A) 

1A9 FORMAT 1/ ,10X .* A3SCLJTE VE_DCITY & ITS DIR. »ITH THE X CJO. '*//) 

150 r0RMAT(//,20X, *1 » .6X.* J ‘.SX.'VV • , 20X . • THEY A* , / J 

200 F0KMATI7X.20X,* I '.SX.'J •.SX.'VVBX • ,2CX. • THETA* .✓> 

300 FQRMAT{//’.20X.*I '.SX.'J •.SX.'VVav • . 20X. 'THET A ' . /I 

RETURN 
END 
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SU3K0JTI HE SHA3AN 
IMPLICIT SEALES (A-H,0-2) 

COMMON /ONE/ VI 

common /t»o/ I lo , i i*r,isym,ipap,icod t ip. XCOM 

COMMON /=OUR/ Rl.RF.lPREVO 

common /dims/ R. ftC.M* e.DX, Or , 81. B6 

common /=0RK/ Ml .M2, IC.OC > I IC. JJC.H3 t N3 

common /MSH/ si. S 2.S3. SA.Pi ,P2,P3,P4 . JU. JU.il. ip 

COMMON /PrtI/ P (pD, 30 J .PSY( 2,30 I, PBX{ 50.21, PP(50, 30}, AS yiI2, 30) , 
1 PBXl £ 50,2) ,Pt 50,3 0) 

common /V£E5/Ut5 0, 30) ,V( 50,30),USX( 50,2) .V8X£ 5C.2J . 

I UBY(2.30}.VBy(2.33>.Ul(50.30).VO( 3 0,30), USX I (50.2) 

I ,vaxl(S0,2) ,UBY1£2,30) ,VBYl£2.3a) . VV(50,30) , ANGLE t 5 0, 30 ) 
CCMMON /CAL/ M2l , ICI .JCl.Ml 1.N31.IC2 
common /3n/ rg .gaka.to.wv ,ve 

COMMON /RO/ ROAt 50,30) .R03X(50.2) ,ROBY(2,30) 

CCMMON /GO/ yVEX(50.2).VVBY(2,30).ANGLl( 50,2 ), ANGl2(2.30) 

CC CALCULATE =OISSONS SOURCE 01 S. Ftl.J) 

C 

DC AOO 1=1 ,H3 
IFdPG.SO.l) GO TO 101 
DO 400 J=1,M2 

400 F(1.J)=O.ODO 

IF( tSYM.EO. 1) IIC=M1 

IFtlSYH.EO. 1 .AND.IPAP.EQ, 1 ) GO TO 777 

AA=O,0D0 

DC SSS 1=2, ICI 

CALL INTPL(I) 

555 AA=AA+( JU-JL)*DY»DX 
DO 5556 I=IC2, Mil 
CALL INTPL(I) 

55S6 AA=AA«-( JU-JL)*DY*DX 
DO 666 1=M1,M3 
CAUL IMPL(I) 

666 AA=AA+{ JU-JL)*OY*DX 

A 1=£M3-1 I C)*CN 3-2 )*OX*DY 
IFdSYM.EQ.l ) Al = (M3-llC)*tN3-l )=40X*DY 
HM=3/( AA-Al ) 

IFUFREVO.EQ.t ) G3 TO 3336 
GO TO 1 
3336 continue 

RI l=Rl +H+R 
DC 401 1=2, Ml 

R22=R1+H+{MI»I )*OX 
CALL 1NIP_(I) 

DO 401 J=JL,JU 

401 F( 1. J)=R11*RH/R22 
GO TO 3337 

777 CONTINUE 

00 11 1=2, Mil 
CALL INTPL(I) 

DO 12 J=JL,JU 
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X=(I-1C)*0X 

V=(J-JC)Oy 

39090 03 »3»( l.ODO-l X»X+y*=rj/{ K4RJ ) 

12 CONTINUE 

11 continue 

3337 Ir {ICOO.EQ.O ) GO TO 7 

1 IFtiCOQ.EO.Ol GO TO 17 
CALL HOSNy 

DO 2 1=2, M3 
CALL INT^LtlJ 
DO 2 J=JL,JU 

2 Pt I, J)=ROA( I, J ) 

DO 3 1=2, M3 

?6X( 1,1 )=ROaX{ 1, 1) 

5 PBX( 1 ,2)=R0BX( 1,21 
DO A J=1,H21 
PBY<l,J)=RDBy( 1, J) 

A P6YI2, J)=R08y<2, J) 

IL0=1 

CALL HVVE- 

IFIIPAP.EO. I 1 R3=0.DD0 

IFUFREVO.EQ.il Rrt=O.ODO 
IFUSyw.EQ.l 1 IIC^Ml-l 
DO 5 I =2, M3 
lFtI.G7.IIC) RH=0, ODO 
CALL 1NT»_<I) 

DO 5 J=JL,JU 

5 Ft J, J)=F< 1 , J)*(ftO3y<2,2)/R0A( 1 , J 1 1 FR h* I ROB Yt 2 . 2 1 /ROA ( 1 ,JJ 1- 
l <U1 { 1 ,J) AUt 1, J) +V0 < 1,J )»V< I . J J )/TOA( 1,J 1 
GO TO 7 

17 DO 6 1=2, M3 

IFtl.GT.IlC) RH=0.000 
CALL INTSLtl) 

00 8 J=JL,JU 
a F(I,J)=R^ 

7 AE=0.000 

DO AAA 1=2, Id 
CALL INTPL(I) 

J1=JL 

J2=JU-1 

DO AAA J=J1 . J2 

AAA AE=AE + ( (=tI,J}+rU + l , J)+Ft 1 , J+ll +Ft I + l ,J+l ll/A 1 + OX ADY 
00 AAA6 1=1C2.M1 1 
CALL INTPLtl) 

J1=JL 

J2=JU-1 

DO AAA6 J=J1 , J2 

AAA6 AE=AE+( ( =tl , Jl+FU-1 , J J +r ( I . J + l 1 +Ft 1-1 ,J+1 1 1/4 lADXADY 
00 333 1=M1,H3 
CALL INTPLtl) 

J1=JL 

J2=JU-1 

DO 333 J=J1.J2 
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333 AE=A£ + ( (FU , J)+F{I-1 . J ) + F 1 1 . J + 1 ) +F( 1- 1 , J + DJ 7A ) :fUX*j r 
AN=(3-N31 + 0Y)*‘V1*RCI5Y{2.2 J 

00 29 J=l.')31 

29 AN=AN+VI*OY«:fiOBYC2.2) 

1 F 1 1 SYM.5'3. 1 .AND. 1 PAP .£0. 1 I AA=1 .5 71 4 3 000*K*!< 

OO 30 1=2. M3 

CALL INT=LII) 

DO 30 J=JL.JU 

30 F t I. J )=F( I . J }-( ah- an l^-AA 
101 CONTINUE 

IFU^R.cD.ll GC TU 300 
»RITE{6 .200 } 

DO 9 Is=2iM3 

9 l^RITEC 6. 100 ) 1 ,( FC I • J) >J=1 .M2U 
WRITEIo . 14) 

1 4 ~CRMAT <///'« ' 

1 »«:» vti**9***** *^*****^*’i‘*** ******* ^* *.*************** • .// J 

100 FCPMAT(//»7X,* X = • . I 10./ . 4 ( 8F10 .5 1 ) 

200 *ORMATC/. 1 OX.* POISSONS SOURCE D1 S, AT EACH GRI O POINT*.//) 

300 CONTINUE 
RETURN 
END 
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SUBROjriNE H3$vr 
IMPElCiT FitAL^S {A-H,0-2) 

CUMHQN /TWO/ IU3.IWft.lSyM,I?AP,IC03.IP.IC0M 

COMMON /FOUR/ Rl.^F.IFREVO 

CCMMOM /DIMS/ R> RC<H. B.DX.DY tSI tBB 

COMMON /FO-iM/ Ml ,M2> 1C < JC •! IC. JJC.M3 .N3 

COMMON /MSH/ SI* S2.S3.SAt?l tP2»P3.PA, JC.JU.Il-. 1ft 

COMMON /VEES/U (50>3U).Vl50f30}>UBX(bO*2I*VBX(50*2]> 

1 UdYt 2.30) ,V3 V< 2*30) » UK 50,30 J , VOISO.30) ,U8X1 <S0,2) 

1 > VSXl ( 53*2) .UBY 1 {2.3C }* V BY 1 (2>30) * VV( 5C.30 ) . ANSl£( 50.30) 

CCMMON /CAL/ M2) *ICl .JCl .Ml 1 *N31 ,IC2 
COMMON /SH/ RE .GAMA.TO.WV, V£ 

COMMON /RO/ RQA(50.3C ) ,R0BX<50.2),a0ayt2,30) 

COMMON /OD/ VVEXt 50.2 ) • VVBYC2 .30 ) . ANGLl ISO .2) , AAGL2( 2> 30 ) 

C 

CC CALCULATE DENSITY ftOAU.J) AT EACH GRID POIN 

CP=RG*GAMA/(GAHA-1.000) 

POY/ER=1.3 00/(3AMA-1 .ODC) 

83 0K=VE»V=/(2.0D0eCP*T0J 

IF (IPAP.EO. 1 ) GO T3 A 
00 1 t=2.M3 
CALL INT^Ld) 

00 I J=J_,JU 

VEL=VV t I .3) *=VV{I , J ) + WV»WV 
T=l. 000-DK*VEL 

IFIT.LE.C.OOCO IDDD ) GO TO 30 
ROAII ,3)=T**PDViER 

1 CONTINUE 
GO TO 8 

A CONTINUE 

DO 5 1=2, M3 
CALL INTPLII) 

00 5 J=JL.JU 
X={ t-IC)*DX 
Y=(J-JC)*DY 

WV=S.09OD0«B«=( 1.0D0/RF)*{1 . QD0-< X*X+ YAY )/ lft*R ) ) 

V£L=VV{ I .J) *VVU . J ) + WV»*V 
T=1 .f)DO-OK*VEL 

IF(T.LE.O.OOOOIODO) GO TO SO 
ROAII, J )=T**POWER 
5 CONTINUE 
8 DO 2 1=2, M3 
DO 2 3=1.2 

V£L=VVBX <I,3)*Wax(I,3) 

T=l.ODO-DK* VEL 

IF(T .LE. 0.00001000 ) GO TO BO 

2 RQBX (1 . J)=T^=fcPOWEft 
DC 3 3=2. M21 

OD 3 1=1.2 

V£L=VVBY(I ,3)AVV3YI1 .3) 

T=1.000-DK*VEL 

IFIT.LE.0.00CQ13O3 ) GO TO SO 
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3 ROav II, J)=T**BOWE5 
GO TO 90 
oO CONTINUE 

*PIT£{o,32> I.J.T.VEL.DK 
S2 FOKMAT(5X,2I10 ,3FI 2.7> 

V£=.90DO*VE 
WRITE(6,3ii) VE 

64 =ORMATC 20X. * VE=' ,F 12.7.//) 

GO TO 83 
90 RETURN 
END 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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SUBROUTINE INTPL U> 

C 

CC SUBROUTINE INTPLCI) PGR DETERHXNiNG BOUNDARY POINTS FOR. A GIVEN 
CC COLUMN l.tl.JLU)} 15 THE lO*ER MOST INTERIOR POINT AND 
CO (I.JU(I)) IS THE UPPER hOST INTERIOR POINT 
C 

IMPLICIT REALMS (A-h.O-Z) 

COMMON /IrfO/ lL0tIhk,I5YH,IPAP,lC0D. IP.iCOM 
COMMON /DIMS/ ft . RC . H . B i DX, DY, B 1 .BB 
COMMON /FORM/ Ml ,M2. I C. JC. I 1 C . J J C. M3. N3 
COKMON /MSH/ S 1.S2. S3. SA, PI ,P2,P3 , PA .JL.JU.il. I ft 
C 

IFU .GT. 1C) GO TO 10 
X=II-1)*DX-R 
YY2=DABSIRAR-X *X ) 

YV=DS0RTIYY2> 

IFlISYH.ED.l) GO TO 9 

T1=R-YY 

Y2=R+YY 

JL=< Yl/DY + O.OlOOODOI+2 
JU=(Y2/Oy-0 .OlOOODOJ+l 
GO TO 40 

9 JL=1 

JU=IYT/OY-0 .01000)41 
GO TO 40 

10 continue 

IF (I .EQ.Hl ) GO TO 19 
IFU.GT.Ml) GO TO 20 
X=U-l J4DX-R 
YY2=DA8S(R*R-X*X ) 

VV=OSORTIYY'2) 

IFdSTM.EQ. 1 } GO TO 18 
JL=2 
Y 2=S +YY 

JU=IY2/DY-0 .OIODQDO )+l 
GO TO 40 
IS JL=1 

JU=( YY/DY-O-Ol ODO)+1 
IFIYY.lt. B) JU=N3 
GO TO 40 
20 CONTINUE 

IFUSYM.EQ.t) GO TO 30 
JL=2 

IFII.GT.IIC) GO T3 30 
X=RC-( I -Ml )*DX 
y Y2=0A3S( RC*RC-X*X) 
rY=DSQRTtYY2) 
ri-( RC+B)-YY 
JU=l Y1/DY-0.01000DO)+1 
GO TO 40 
1'9 CONTINUE 

IFt tSYM.EO. 1 > GO TO 30 
JL=2 


1 

2 


5 

o 


9 


10 

11 

12 

14 


13 

17 

19 


20 

21 

24 


27 

29 

30 

31 
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JU=JJC 31 M 

GO TO 40 32 

30 CONTINUE 53 

JL=2 3t 

IF II SYM.EQ. I ) JL=1 
JU=l B/OY-0 .001 OODOlt 1 
4 0 CONTINUE 

RETURN 
ENO 


ORIGINAL PAGE IS 
OP POOR QUALITY 
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SUBROUTINE MEaHU.J) 
implicit ^EAL*8 (A-rt.D-Z) 

COMMON /mo/ ILO. IWRt ISYM, IPAP, ICHO. IP. ICUH 
COMMON /OIMS/ R.RC.h.B.OX.OY.ai ,63 
COMMON /FCRM/ MI,M2,1C.JC.IIC,JJC.M3,N3 
CCMMON /.MSH/ S l,S2,S3,SA.Pl.P2.P3. = A, JL.JU.IL.IR 

COMMON /o-U/ PIS0.3G) ,PBY(2.30J .PBX150.2) ,PP I 50 .3u J , PBY 1 12. 3u). 

I PBXI <50,2), Ft 5O.30J 

C determine the standard five points LAPLACE DIFF. OPERATOR AND 
CC THEIR POSITIONS S1.S2,S3,SA a.R.T, POINT O(X.Y) 

C 

IF( J.EO.JL.ANO.J.EO. JU) GO TO 101 


IFIJ.EO.JL) GO TO 50 COl 

IFIJ.EQ.JJ) GO TO 53 002 

52=1.000 

P2=?<I,J+1) ^ 004 

SA=l.ODO 

PA=?{1.J-1) QOo 

GO TO 70 C07 

SC continue 00b 

IFtI SYM.EO. I ) GO TO 51 

IFU.GT.IC) CD TO 51 009 

S2=l .ODO 

P2=PU,J+1) 0011 

X = R-U-1)*DX 0012 

yY2=DABS<R*H-X4X) 

YY=DSQRT( YY2) 

Yl=YY-( JC-J)*DY 
S4=Y1/0Y 


IFISA.GT.l.ODO) SA=1.0O0 
P4=P3X ( I .1 ) 

GO TO 70 


51 CONTINUE 0019 

S2=l .300 

P2=P{I,jH-n 0021 

SA=1 .000 

P4=pax£l,l) 0023 

30 TO 70 0024 

60 continue 0025 

IF(I.EQ.Ml) GO TO 52 0037 

IF(l.GT.Ml) GO TO 65 003fa 

X=CIC-I )*OX 
YY2=DABSCR*R-X*X) 
yY=osoRT<yy2) 

Y1=YY-£J-JC)*DY CC30 

SZ=Y1/DY 0031 

IF£YY,LT.B) S2=£ 3- ( J- 1 ) *D V ) /DY 
IF£S2.GT.l.0D01 S2=l .ODO 

P2=PaX(I,2) 0032 

S4=l .000 

P4=P£I,J-1) 0034 

SO TO 70 0035 
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62 CCiNTInUi 

I F( ISVM .=3 . 1 ) GO TO 64 
S2=l .DD3 
P2=f>3Xt I »2) 

S4=l .000 
P4=P( I , J-1 ) 

GO TO 70 
6£ CONTINUE 

IF ( ISYH.fcO, 1 ) GO TO o4 
IFU .GT.I IC) G C TO 64 
XC=RC-t I-Ml )40X 
YY2=UA3S{RC4RC-XC*XC I 
YY=DSORT( YY2J 
YC=B+RC-YY- ( J- U»OY 
S2=YC/0Y 

IF CS2.GT. 1. ODO 1 62=1.000 
P2=PBX{ 1.2) 

S4=l .000 
P4=PU ,J-1J 
GO TO 70 
04 CONTINUE 

Y2=9-{ J-1)*DY 
S2=Y2/0Y 

IF ( S2.6T.1 . ODO } ^2 = 1 .000 
32==BX( 1,2) 

S4=I .ODO 
P4=PCI . J-I ) 

CO TO 70 
101 CONTINUE 

Y=6-( J-l ) 40 Y 
5 2=Y/DY 

lFtS2.GT.1.0DD > S2=l .ODO 
P2=PBX{ 1,2) 

S4=l .000 
P4=P3X{ I . 1 ) 

70 CONTINUE 

IFU .Ea.It.) GO TO 80 
IF(I.ECi.lR) GO TO 90 
61=1.000 
PI=P(I-1,J) 

63=1 .000 
P3=P(X+1 . JJ 
GO TO 100 
80 CONTINUE 

Y=( JC-J)*OY 
XX2=DABSt R*R-Y»Y) 
XX=DSQHrtXX2) 

Xl=XX-t IC-I >4DX 
5I=Xl/DX 

IFISl .GT.I. 0001 61=1.000 
Pl=P3Y( l.J) 

53=1 .000 
P3=P(I+1 . J) 


oeiowao k 

OS' JOOR aUAim 


0048 


6060 

0052 

6063 

0054 

0u56 

0056 


0059 

0050 

COol 

0063 

00»4 

006S 

0065 

0067 


00o8 

0070 


0071 

0073 

0074 

0076 

0076 

0079 

0060 


ooss 

0086 

0087 


PAoE 0002 
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GC TO 1 00 
SC CONTINUE 

lf«J.GT.N3) GO TO 91 
Sl=l ,0D0 

Pl=3( I_l, J) 

XX=H-tM3-MI J*OX 
S3=XX/DX 

XP<S3 .GT . 1 . ODO > 53=1.000 
P3=PBy (2, Jj 
GO TO 100 

91 CONTINUE 

IF C ISYM.EQ. 1 ) GO TO 93 
IFlJ.GT.JJC) GC TO 92 
S1=I .000 
PI=P(1-1 . J) 

YC=3*-KC-( J- 1 }• DY 
XXC=OABS ( RC ARC-YCFYC I 
XC=OSaRT(XXCl 
X1=RC-XC-C 1-Ml )»OX 
S5=X1XDX 

IF (53.GT .1 .000 J S3=1.0DO 
P3=PBY( 2, J» 

GO TO ZOO 

92 CONTINUE 
IFIJ.GT.JCJ GO TO 93 
Sl=l .ODO 

P1=P(I-1 ,J) 

S3=l .000 
P3=PSY (2. JJ 
GO TO 100 
93 CUNT INUE 
Sl=l .ODO 
P1=PU-1,J) 

Y={J-JC)*OY 
XX2=0ASS(R*«-Y4>Y ) 
XX=0SQRT{XX2J 
Xl=XX-t 1-1C)*DX 
S3=X1XDX 

1 F < S3.GT. 1 . 000 ) 53=1.000 
P3=P5Y(2.J J 
100 CONTINUE 
RETURN 
END 


0090 

00101 

00102 

00104 

00109 
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00116 

00117 
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subroutine IRLCJ) 

IMPLICIT BEAL*S lA-H.O-Z) 

COMMON /T*OZ IU3»1 WR , ISYM. IPAP. I CUO , iP.ICOM 
COMMON /DIMS/ R.BC,ri.B.OX,Dy,BI,53 
COMMON /form/ hi ,M2 , IC.JC, I IC» JJC» H3 . N3 
CUHMON /MSN/ Sl,S2tS3.SA.Pl.P2tP3.PA« JL.JU.JC.IR 
C 

CC SUBRCUTINc IRLIJJ FDR OETEPMIMNG BOUNDARY POINTS F05 A IVtN 
CC ROM J,(J.I^(J)) IS TNE FIRST INTcRIQR POINT tiJ.IKIJII IS THE MAXIMUM 


CC INTERIOR POINT 
C 

IFCJ.G7.N3) GO TO lOOO 01 

Y= (JC-J )*DY 

XX2=DABS(RT=H-Y*Y) 

XX=DSORT( XX2J 

X1=K-XX 05 

IL=( Xl/DXFO.Ol OOQOOI+2 

iR=H3 07 

GO TO ACOO Ofc 

ICOO CONTINUE 09 

IF(ISYM.E0«1) GO TO 3000 

1F(J. GT.JJCJ GO TO 2000 010 

Y=R-{J-I)AOY Oil 

XX2=DA3S{RPR-Y*Y) 

XX=DSQRT (XX2I 

X1=R-XX 014 

IL=I X173X+0 .01 00000)42 

YC=a+RC-{ J-1 )*0Y 016 

XXC=DA8S{HC*RC-YC* YC) 

XC=OSQRHXXC) 

X1=RC-XC 01 St 

IP=( Xl/OX-0.01 00030)4N1 

GO TO 4000 021 

2000 CONTINUE 022 

IF(J.GT.JC) GO TO 3000 023 

Y=< J-1 )*DY-R 024 

XX2=DA3S{RT=R-Y4YI 
XX=DSaRT{XX2) 

X1=R-XX 027 

1U=( X17DX+0. 0100 000)42 

IR=«l-l 029 

GO TO 4050 030 

3000 CONTINUE 031 

Y=r{J-JC)t*DY 
XX 2 =DAeS{R 4 fi-Y*Y) 

XX=950RTtXX2) 

X l=R-XX 035 

IL={ X1/'DX4 0 .0100 000) 42 
. IR=tXX/DX-0.01 OOOOOJ41C 

4000 CONTINUE 040 

RETURN 
END 
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h.UaROUTINE ECSIdl 
IMPLICIT ^£AL*6 (A-rl.O-Z) 

COMMON ^TtlO/' ILD . I*R« ISTM, IPAPi ICOO* IPt ICOM 
COMMON /DIMS/ R.8C tOy ,ai .83 
COMMON /FORM/ Ml ,M2,1C.JC.1IC. JJC.M3.N3 
COMMON /MSri/ S J, S2,S3,oA. =I .P2.P3.=4, JL, JU. IL. 18 

common /P-U/ P t Sli.30) .PSVt 2 .30) .P8< (SO .21 .PPI&0.30), P3V I 12, 30 J , 
1 P3X1 (SO, 2>.F< 50,30) 

COMMON /INPC/ PL.PU.pt. PR 


CC CALCULATE POTENTIAL =UNCT1DN AL3NG THE B. FOR EACH X {P3XI.P8X2) 
C 

IF(1 .68. XC) 60 TO 180 
J=JL 

CALL IR_(JJ 
CALL MESH(I.J) 

X=(IC-l )*3X 
yy=DABS<RPR-X*X) 

Y=0SQRT( yy) 

1F< ISYM.Ea. 1 )' GO TO 100 
S5=SA»Dy*X + 0X(‘Y 

PL={ 1+SA)TP( 1+ I, J)-5A»P{ I+l .J+1 ) 

PSXI I. I )={ SA+? Y*X*?_+DX»y*P( X , J) )/SS 
GO TO 101 

100 P3X( I . I)=PU .J) 

PL=P3X (1,1) 

101 J=JU 

CALL IRL(J) 

CALL MESH(I.J) 

SS=S2*^DY*X+DX*Y 

PU=( 1+S21 *P( 1+ 1. J)-S2*PCI 4-1 , J-1 ) 

PSX( I. 2)=( OX»y*P { I ,J)+DY*S2*XTPU)/SS 
GO TO 900 
ISO CONTINUE 

IFU .GT. ICJ GO TO 190 
J=JL 

PSX( I , I ) = P( I . J > 

PL=PBX (1,1) 

J=JU 

P3X( 1, 2)=P{ I .J ) 

PU=PBX(1,2) 

GO TO 900 
190 CONTINUE 

IFU .Gc. Ml ) GO TO 210 
J=JL 

PBX< I, 1)=P( I.J) 

PL=PBX ( I . 1 ) 

J=JU 

IP1J.LE.N3) GO TO 218 
CALL IRL(J) 

CALL MESHd.J) 

X=( I-1C)«X 
YY=DABS (R*R— X#X) 


1 

2 

3 

A 


7 


11 

12 

13 

^14- 
1 5 
lo 
17 
lb 

19 

2C 

21 

22 

23 

2A 

28 

2b 


27 

28 


41 



s IV G 


_EVEL 21 SCSI 

y=OSQRT lYV ) 

SS=D Y»i>2*X+DX* Y 

PU=t 1+S2)=*‘P( I-l, J)-S2»PU - 1 . J-l) 
PBXd t2) = {DX«Y VPU , J)+S2*0Y+X*=PU)/i>S 
30 TO 900 

21o PBX( I , 2 )=“( I , J ) 

PU=P3XtI ,2) 

GO TO 900 
2X0 CONTINUE 

IFdSYM.EO. 1 ) GO TO 214 
XF(1 .GT.Ml) GO TO 212 
J = JL 

PEXIl ,1)=P(I,J ) 
pi.=paxi 1 , 1 j 
j=jj 

PBXd ,2)=Pt X-1 .J+1 ) 
pu=paxd . 2 } 

GO TO 900 
212 CONTINUE 

IFU.GT.dC) GO TO 214 
J=JL 

PBXdil J=Pd.J J 
PL=PBxi I . n 
J=JU 

CAI.L IRLIJ) 

CALL MESHd,JJ 
X=RC-d-Ml J *DX 
YY=£) ASS t RC*RC-X4X> 

Y=DSORT( YY) 

SS=DY*S2*X+0XF Y 

PU=d + S2)*P{ I-l , J)-S24Pd-l , J-1 ) 

^bXI X.2) = OY4S24X*PU+DX4Y4Pd • J) ) /’SS 
GO TO 900 
214 CONTINUE 
J=iJL 

PBXd, 1 ) = = ( 1. J ) 

PL=PBXd ,1 ) 

J=JU 

P8Xd,2) = = ( I. J J 
PU=PBX( 1,21 
900 CONTINUE 
RETURN 
End 


DATE = 772e>J 21X45/4 1 PAGE 


31 

32 

33 

34 


35 

3o 

37 

35 

39 

40 

41 

42 

43 

44 

45 

4d 

47 

4b 


bX 

52 

53 

54 


5b 

57 

Sb 

59 


0002 
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SUoRODriNc BC3J{J>V1> 

1^*PLIC^T HEAL»e (A-H,0-2) 

COMMON /T«0/ ILO. ISYM, I=AP,XCQDtlP»ICDM 
COMMON /DIMS/ B. MC.H.B.OX.OY.Bl.aa 
COMMON /=ORH/ Ml ,M2. ICiJC.IlC, JJC.M3iN3 
COMMON /MSH/ 5 1. S2>S3>S4.^ I *P2>P3. 34. JL. JU> { L. IB 

COMMON /PHI / P(5 0,33 ) ,PBV (2 ,30) ,P3X I 50,2),PP150,30J.3 3Y1 12.30 ) , 

1 PBXII50.2), FI 50,30) 

CCMMON /INPO/ PL.PU.PT.PR 


CC CALCULATE POTENTIAL FUNCTION AwDNG THE B. 
C 

IFIJ .GE. JC) GO TO 200 
I=IL 

CALL INTPl(I) 

CALL MESHtl.J) 

Y»I JC-J)*3Y 
XX=DABS ( R*R-Y*Y) 

X=DSQRT(XX) 

SS=OY*X+DX*Sl*Y 

PT={ 1+Sl )»P( I. JTl )-S IPPt I+l. J+ 1) 

PBYl 1,J>=(0 YTX^PII ,J)+DX»SH=Y*PT)/SS 
1 = IR 

CALL INTPLII) 

CALL HESHll ,J) 

IF(u .GT . N3) G3 T3 110 
PSYI2, J)=Pt I , J )+DXP33*Vl 
PR=P3Yt 2, J) 

GO TO 900 
1 1 0 CONTINUE 

IFIJ .GT. JJC) 30 TO 120 
Y=e+RC-(J-l J *DY 
XX=OAaS< RCPRC-Y^Y) 

X=DSQRTI XX ) 

SS=0X*S3*Y+DYTX 

P»=(l+S3)*PtI . J-l)-S3TP(I-l ,J-1) 

PBYI 2. J)=I D¥*X*PI I ,J)+DXTS34>Y*PR)/SS 
GO TO 900 
120 CONTINUE 

PSYC2.J)=3(I,J1 
PR=P8Y (2, J) 

GO TO 900 
200 CONTINUE 

IFIJ .GT. JCJ GO TO 210 
1=IL 

P8YI I, J)=3£ I . J ) 

PT^PBYI 1 . J) 

I=IR 

CALL INTPLII) 

CALL MESHII . J) 

IFIl SYH.Sa. 1) GO TO 101 
PBYl 2.J)=3{ I.J ) 

GO TO 102 


FOR EACH J IPBY1.PBY2) 
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101 PBY(2-.JJ=P( I . J )+D<»33*Vl 

102 P«=PdY(2.J) 

SO TO 900 

210 cgntinue 
I = U- 

CALL INTP.d) 

CALL MESH ( I t J) 

Y=( J-JC) *DY 
XX=DAaS ( R*R-YEY ) 

X=OSQRT(XX) 

SS=DX*S1*Y+D Y*X 

PT=( 1+Sl ) tP( I . J- ll-51*PCI+l. J-lJ 
PSY( 1. J3=<DX*Sl*y*PT + DY*X*P(l. Jj >/3S 
r=iR 

CALL iNTPLd ) 

CALL MESHd .J) 

IFCISYM.EQ.03 GO TO 260 
IFCJ»3T.N33 GC TO 260 
P0YC2.J3=^C1 »J 3+DXTS3* VI 

PR=PaYC2. J) 

GO TO 900 
260 Y=IJ-JC3I'0Y 

XX=0 A3S ( YT^Y) 

X=DSQRT{XX) 

SS=DX«S3*^Y+DY* X 

PP=( 1+S3 3*P( I . J-1 J-S3*PC J- 1. J-1) 
PBYC2. J)=( OY*X*P(a , J)+0X*33»Y*PK)/3S 
900 CONTINUE 
RETURN 
END 
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DATE 
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2 1VAS-/AT 


PAic UUOI 


3U3HQUriN£ riVVEU 
implicit :iEAL»8 (A-HtO-Zi 
COMMON /ONE/ VI 

COMMON /T*0/ ILU.JhR.ISYM.IPAP.ICOO, IP. ICGM 
COMMON /OIMS/ R.RC.-j .B.DXjOY.BI.OB 
COMMON /FORM/ Ml . H2. I C. J C. 1 1 C. J JC. M3 . N3 
common /MSrI/ S 1.32. S3 ,SA .PI .P2.P3. PA .JL.JU.il. IN 

COMMON /PHI/ =(50.30) .P3y( 2.30) .P8X( 50.2) .PP(50.30) .P3Y1 (2.30) • 

I PBXl CS0.2).F(S3.30 ) 

COMMON /VEeS/U(50.30) . V(5 0 .3 0 ) ,U8X( 50 .21 .VBXISO .2) > 

1 UBYI2.30). V3Y( 2.30 ) ,U1( 50.30) . V0{ 50,30 ,USX1 { SC. 2) 

X . Vaxi (50.2) ,UBYI (2 ,30) .VBY 1 ( 2 ,30 ) . V V ( 50 ,30 ) . ANt.LE( 50,30) 

COMMON /INPO/ PL.PU.PT.PR 
C 

CC CALCULATE VELOCITY COMPONENTS J.V IN THE X£Y DIK. (U=0=/DX ,V=DP/JY) 
C 

IF( Iiim.EQ. 1 ) GO T3 25 
WRITE(6,14) 

WRITE(6.400) 

WRITE(6. 107) 

25 CONTINUE 

DO 10 1=2. M3 
CALL INTP_(l) 

DO 20 J=JL.JU 
CALL IRL(J) 

CALL MESH(I.J) 

C13=51/CS3+(S3+Sl) ) 

C31=S3/{S1*< S3+S1) ) 

C1=<S1-S3)/(S1^S3) 

C24=S2/(SA*iS2+SA) ) 

CA2=S4/( S2*( S2+SA) ) 

C2= ( SA-S 2 ) / ( S2 *5 A ) 

UCI, J) = {P3*C13-Pl*=C3l-P(l . J)AC1)/DX 
IF(J.EQ.l) PA=P2 

V( I. J)=(P2ACA2-PA=C2A-P(I, J )*C2)/0Y 
IF(lLO.EQ.l) GC 70 20 
01(I,J)=J{ I.J) 
ve (I ,J)=V(I . J) 

IF(IWR.EQ.l) GO TO 20 
kR1T£(6.111) I .J.J(1>J),V(1,J) 

20 CONTINUE 
10 CONTINUE 

IF( Xi^R.EO. I > 30 TD 30 
IFdoYM.EQ. 1 ) GO TO 656 
WR1TE(6, lOS) 

GO TO 30 

566 k.RITE{6.1 09) 

30 IFUSYM.EO.l) I1C = M1 
DO 50 1=2, M3 
CALL INTPl(I) 

IF(1.EQ.M3) go TO 100 
J = JL 

CALL IRL(J) 
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CALL MES^( I , J) 
caul 3CSI [ I ) 

IFU.GE.ICI = i + l , 1 ) 

UBX( I .1 )=C PL-PBX ( I ,1 I jyox 

vsx( I, I )=( p ( I , j)-P3xi 1 , u J /csA*Dy) 

J=JU 

CALL IRL(J) 

CALL MESHtl.J) 

CALL 0C3 1 C I ) 

IFU.GE.IIC) PU=P3X( I + l ,2) 

JBXI I. 2 >=l PU-PBX(I ,2)) /OX 
IFU.3T.ICJ UBX( I,2I = -UBXU,2) 

JFa.GE.ltC) USXd ,2)=-UBXU .2) 

Vax(l.2> = CPBX{l,2)-P(I,J))/{ S2»DY) 

CL TO 900 
100 CONTINUE 
J=JL 

UBXC 1 « 1 )={PSX£ I. 1 )-PBX(I-l , 1 J )/OX 
VBXt t . n = tPt I 1 JJ -PBXt 1,1) )/DY 
J=JO 

CALL IRLtJ) 

CALL MESHtI»J) 

CALL BCSKl ) 

USXU ,2) = tPax( 1,2)-PBX(1-1,2J J/DX 
V5X(1 ,2)=[PBX( 1.2J-PII ,J))/tS2*OY) 

900 continue 

IFdLO.EQ.l) GO TO 50 
uaxui , l)=UBXl 1, U 
VBXl 11,1 ) = VSXt I, 1 ) 
uaxi ( 1 .2 >=uex( i ,z> 

VSXK I , 2)=VBX{ I , 2) 

IFUlvR.EO.I) GO T3 SO 
IF(1 SYM.EQ. 1 J GO IQ 103 

WKITEIB, 310 ) UBX( l.l),VbXU,l> .UBXt 1,2) .VBXIJ ,2) 

GO TO 50 

103 WR1TE(6,210) USXtl ,2 ) , VBXl I , 2 ) 

50 CONTINUE 
M21=M2-1 

' IFU <(R. =0. I ) GO TO 15 
UR ITEt 6,777 ) 

15 CONTINUE 

00 oO J=2,H21 
CALL IH_tJ) 

I=IL 

CALL I NTPut I ) 

CALL MESHCI,J) 

CALL BCSJt J, VI ) 

0BY{ 1, J)=i P I I , J)-PBY( I , J) ) /< Sl*DX) 

VBY( 1, J )={PT-PBYt 1, J ) )/0Y 
IFJJ.GE.JC) VBYCl ,J )=-VBY( 1 . J ) 

1 = IR 

call INTPL;Ci1), 

CALL MESrt(I.ij) ' 

. ' • t - ! 
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CALL BCSJI J.Vl ) 

IF(J.ST.N3) GO TO laO 
USYC2. J > = (PBY{ 2, J)-Pt I. J ) )/( S34DXJ 
V3 Y{2. J)=0. OOO 
30 TO 3000 
280 CONTINUE 

IFdSYM.EQ. 1 1 GO TO 240 

IFIJ.GT.JJCJ GO TO 190 

UBY(2,J1 = {=BY(2. J)-P( I.JJ)/{i3FDX) 

VBV(2. J>=(PSY<2t JJ-Pftl/DY 
GO TO 2000 
190 continue 

IF(J.GE.JC) GO TO 220 
UBY( 2. J) = { 3SY( 2. J)-P( I ,J))/DX 
VBYta, J) = (PBY(2, J+1 )-PBY(2»J J )/DY 
GO TO 2000 
220 CONTINUE 

XF(J ,GT. JC) GO TO 240 
UBYI 2, J) =<PBY( 2. J)-P(l . J1 )/OX 
VSY{2. J)=(P8Y{ 2. J)-PBY(2, J-11 J/DY 
GO TO 2000 
240 CONTINUE 

UBY(2. J}=t33Y(2. J)-P< 1. J ) l/( S34DX> 

VBY{2, J}=(P3Y{2,U)-PR)/DY 
' 2000 CONTINUE 

IF(ILO.EC.I) 60 T3 80 
UBYI (1 ,J)=UBY( 1. J) 

VBYK 1 . J>=VSY{ l.JJ 
UBYI C2.J)=UBYt 2. J) 

VBYl 12 .J»=VBY( 2. J) 

IF(IMR.EO.I) go TO 60 

WRITE{6>313) JBY(l<Jl»V8Y{l.Jl»UBY(2>JJtVBYI2<U} 

60 CONTINUE 
14 FORMAT! IHIJ 

ICT FORMAT!//, lOX, *I •.SX.'J ‘tSX.'U •,20X,*V *J 

1 08 FORMAT!//. 10X,*uaxl •.ISX.'VBXl ',lSX,*UaX2 »»1 SX,‘VdX2 •} 

111 FORMAT! lOX. IS. 5X.1 5, 1 OX. 31 5.6 . 1 0 X. 3 I S . 6) 

210 FORMAT!//, 10X,2(SX. 015.6)1 
310 FORMAT!//, lOX. 4! 5X. 015. 6)1 

400 FORMAT ( lOX. •VELOCITY COMP. U&V IN THE X6 Y DIR.*.//) 

777 FORMAT!//,! OX, *UBYJ •,15X,»VBYI • , ISX, *UBY2 •, 15X, •VBY2* ) 

109 FORMAT! //, I 5X, *UoX2 '.ISX.'VBXE ■,/) 

RETU RN 

END 
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SUBkOUTI ME PLTISCH JMAX C, I S Y H3L, 1 UPOD W, b 1 2E. VISOBR t VAR . X. Y ) 

IMPLICIT RSALsa lA-H.D-Z) 

CCMMON /MSH/ SltS2.53.SA,Pl,?2.P3t^A, J1.,JU,IL, IFi 
31 MENSICN VAR< 32 >21} >M32>21).r(32>21} 

1C=ISYMBL 
I CA=IUPQD»i 
JMAXCl = J MAxC-1 
NNN=1 

1115 CONTlNJE 

1113 3C 1128 J=1.JHAX01 
CALL IRL(J) 

DO 1130 INV=IL>1R 
IIB=1 

I=1R+I_-INV 

1103 IF( (VARd > J) .GT. VI S03R) .AND. ( VAMU . J + 1 ) .GT . Vi SOBR ) } GO TO 1117 
IF t t VAR( 1 , J 1 .LT. VI bOBK) .AND. I VAR(1 , J+1 } .LT . V 1 sCBR J ) GC Tu 1117 
I IS=1 

DS= XU , J+1 >- X( I . J) 

DSS=0S*{ VISOBft-VAktl , J) ) /I VAR{ I . J+1 )-VARU.JJ) 

C YYB= { YtI,J)+ YU.J + UJ/2. 

DN= YU f JU )- YU , J) 

UNN=DN*{ VISDBR-V AKU ,31 }/t VAft( I . J + 1 1 -VAR( I . J) 1 

XXB= { Xd.Jl + DbSl 

YYe= ( V d , J+1 )-ION-DNN) J 

i*FiirEI6,lllA) XXB. YYB, l.J, VISOBR.VARU ,J+l } . VARU > J) .OS.ON 
1 1 1 A FORMAT (2X,"XXB=*,rl3 ,b.2X, 'YYBs’.FlO.B.aX, *1=*, 13.2 a, *J=*,I2,2X, 

1 > VI SD5R=* ,F 14.6 ,3X, A(F1 0.6 J 1 

IFINNN.EQ.il GO TO 1116 
CALL SYM3QLIXXB,YYB,SI2E.IC.0..-I ) 

Z CALL SYMBOL (' YYB, XXB , SI ZE . I C , O . I 1 

CCC CALL SYMB0_ £XXB. YYB, SIZE, IC.0.,-1CA1 
GC TO 1120 

1116 NNN=0 

CALL SYMBOL tXXB, YYB. SIZE, IC .0.,-l) 

C CALL SYMBOL (- YYB.XXB .SIZE , I C, 0 , ,-l J 

GO TO 1120 

1117 IIB=2 

I F( ( VARU , J 1 ,GT. VISQBRl .AND, (VARU-1 , J) .GT.V ISUBR) 1 GO TO 1123 
IFdVARI I,J1 .LT.VISOBFU.ANO. IVARU-1, J] .LT.V1S05R11 GO TO 1123 
11C9 DN=( YU ,J>~ Y U-I ,J 1 1 

DNN=DN*t VI SOBR-VARU-t ,J) 1 / ( V AR U , J } -V Afi U-l . J 1 1 
YYB= I Yl I, Jl-(DN-ONN) 1 

0S= XU , Jl- X( I-l , J) 

OSS=OS» ( VISCdR-VARU-1 ,J) 1 / ( VARU, Jl -VARU-1 , J 1 1 
XX8= ( X d . J 1- t DS-OSSl 1 

MRITEI 6, ! 1141 XXB. YYB.I.JfVISOBR.VARI 1-1.J]>VARU,J1, 05, ON 

IFINNN.EQ.l) GO TO 1118 

ICA=2 

IFCUB.EQ.2) ICA=1 

CALL SYV33LlXXe,rYB,SIZE,IC.O. ,-ICAJ 
C CALL SYMBOL { - YY B. XXB ,S IZE . I C . 0 • 1 1 

IFUia.EQ.il CO TO 1127 
GC TO 1124 

ORIGINAL PAGE ^ 
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Ilia NNN=0 

CALL SYMBQ_ (xxa. rra. sizti 

C CALL SYM30L {- rVB.XX S.S I ZE ,1 C . 0 . .-1 J 

IFIiia.EO.l) GO TO 1127 
GO TQ 112A 

1120 IF 1 < VARtI-1 . J»-l) -.ST-VlSaSR) -AND. I V AR ( t , J + 1 1 . GT -V ISOB? 1 1 GO 1C 
I 1 12A 

IFKVARI I-l . J+IJ -LT.V JS03K) -AND. ( VARt i - J + I J -LT . VI SCSR ] 1 GO 70 
1 1 124 

1121 0N={ Y(I,J + D- YH-l.J+lI) 

ONNsONi (Vise BR -VAR ( I-l.J +1 ) ) / ( VAKU . J +1 )- VARI 1- 1 , J4-1 1 ) 

YYB= ( Y(I .J+1)-IDN-DNNJ ) 

DS= X(I,J+l)- X(I-1»J+1> 

OSSsOSiRI VlSQBR-V AH t I-l.J + 11 )/ tVAR( I. J+1 )-VAR £ I-l, J+1) ) 

XXB=£ I X(i .J-t-I )-{QS-3SSJ ) 

ViRITEt 6, 1114) XXB. YYS»1,J. VISOSR. VARt I-l . J+I j . VAR< I »J + J J .L)S,ON 

IFCNNN.EQ.l) GO TD 1122 

CALL SyMB0L{XXB.YYB.SI2EtIC.0.*-2 ) 

C CALL SYMBOL (- YYB. XXB . SIZS . I C . 0 . ,-l J 

CCC CALL SYMSOLtXXB.YYBfSIZfc.lC.O.i-ICAl 
GO TO 1127 

1122 NNN=0 

CALL symbol tXxejYVB.SI2E.IC.0-.-lJ 
C CALL symbol {- YYB, XXB. size .1C. 0- .-I ) 

GO TO 1127 

1124 IF( ( VARt I-l , JJ .3T.VIS06H) .AND. { V ARt 1 - 1 . J +U - GT -VJSOBR J -ANQ.t I JB.EQ 

1.1) ) GO TO 1 109 

IFt t VARt l-I » J) .LT. VI SOBR) .AND, ( VAR( I-I»J+I).lT . VISOBRJ .ANO.illB.EQ 

1.1) 1 30 T3 1109 

IFt t VARt I-l . J) .GT. VI SOBR) . ANO- ( VARt 1 - I - J+1 ) -GT . VISOBR) -AND. t JIB.EQ 

1.2) ) 30 TO 1121 

IF ( t VARt I-l . J) -LT. Vl SDB k) -AND- { VAR t 1-1 , J+1 ) -LT .VIS GBR] .AND. t 1 IB.EO 

1.2)) GO TO 1121 
GO TO 1125 

1123 1 1 e= S 

IFt (VARt I-l .J) -GT, VISOBR). AND .( VARt I- 1. J+1) .GT. VISOSR) ) GO TO 1127 
IFt ( VARt I-l ,J) .LT. VISOBR) -AND. tVARt I-l, J+1 )._T -VISGSR ) ) GO TG 1127 

1125 DS= Xtl-l.J+D- Xtl-l.J) 

DSS=DS*( VISOBR-V ARC I-l. J D/tVARtl-l.J+I )-VARt I— I . J J) 

XX8= I Xt 1-1 , JJ+D3S) 

0N= Ytt-l,J+I)- Ytl-l.J) 

DNN=DN» t VXSOBR-VARd-l.J) )/t VARt I-l , J + I )-VAR 1 1-1 . J ) ) 

YYB= ( Y( 1-1 . J+1 J-tON-DNN)) 

»RITet6, 1114) XXB. YYB. l.J, VISOBR. VARt 1- I. J) . VAR tl- 1. J+1 ) .OS, ON 

IFCNNN.EQ.l) GO TO 1126 

ICA=2 

IFtllB.EQ.S) ICA=1 

CALL SVM30HXXB.YY3,SIZE,IC.0..-1CA> 

IF(liB.E3,5) GC TD 1121 
GO TO 1127 

1126 NNN=0 

CALL SYM30L ( XX3 -Y YB . S tZE . I C . 0 . .-1 ) 

C CALL SYMBOL (- YYB. XXB . SIZE. IC . 0 . .-1 ) 
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IV c. 


_=VEL 21 PLTISO DATc. = 77263 21746 /'h1 PAuE OOOo 

IFUIQ.EQ.5) GO T3 I12X 
1 1 2 7 C DNT I NOc 
1130 CONTINUi 
1 12S CCNT INUE 
NNN=1 
RETURN 
END 
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INPUT DATA 


In the f ollowing./ , ,the. different- -parameters' 'used to specify 
the cases that can be computed using the computer program, 
will be explained, then a description in details of how to 
prepare the input data will be given. 

Control Parameters for Computations 

1. (a) Symmetric scroll cross section, ISYM = 1. 

{b) Nonsymmetric scroll cross section, ISYN = 0. 

2. The type of through flow profile is specified using the 
parameters IFG, IFREVO, and IPAP. 

(a) IFG = 1 , IPAP = 0 IFREVO = 0 , for an arbitrary source 
distribution. In this case the source strength is fed 
as an input at all the interior mesh points . 

(b) IPAP = 1, IFG = 0, IFREVO = 0, for a circular paraboloid- 
source distribution . 

(c) IFREVO = 1, IPAP = 0, IFG = 0, for free vortex source 
di s t ribut ion . 

(d) IFG = 0 , IPAP = 0 , IFREVO = 0 , represents unifom 
source distribution. 

3. Either compressible or incompressible flow solutions are 
obtained by specifying the value of the parameter ICOM. 
Compressible flow solutions specify ICOM = 1, and for 
incompressible flow case ICOM = 0 . 

I 

How to Prepare the Input 

The input is divided into seven sets , and is given in 
the following. 

First Set (Fluid Properties) , one card 
READ: RG, GAMA, TO, VE 

according to format (4F10.0) 

Second Set (Scroll Geometry). , one card 
READ: Ml, M2, DO, Rl, B, E, R, RC 

according to format (.2110, 6F10.0). 


ORIGINAL PAGE IS 
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Third Set (Control Parameters) , one card 
5T5AD; ICOM, ISYM, IPG, IPAP, IPREVO 
according to format (5110) 

Fourth Set (Nximerical Parameters) , one card 
READ: ITMAX, EPSMAX, WO, WOP 

according to format (IlO, 3F10.0) 

Fifth Set (Output Control Parameters) , one card 
READ: IP, IVPLOT, IPLT 

according to format (3110) 

Sixth Set (Values for Velocity Potential Contour Plotting ) , 
munber of cards is equal to IPLT 
READ : VISOBR 

according to format (FlO . 0 ) 

VISOBR Nxmierical values of the velocity potential 
contours to be plotted as output. 

Seventh Set (Arbitrary Source Distribution) 

Is required only in the case of arbitrary source dis- 
tribution, i.e., IFG = 1. The value of the source 
distribution F(I,J) is read in DO loop according to 
format (8F10.0) . The input data is fed starting from 
I = 1 to I = M^/ and marching in J direction from 
J = JL, JU as shown in Pig. 3a, (Ref. 1) 

OUTPUT 

The program output includes a printout of the pertinent flow 
properties at all the grid points every IP iteration and of 
two figures, one for the desired velocity potential contours as 
specified by IPLT and VISOBR and the second showing the velocity 
direction in the cross-sectional plane. 

Samples of the program output are included in Ref. 1. ■ 
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